1 



X 



g factor in a light two body atomic system: 
a determination of fundamental constants 
to test QED 
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j^y^ ' Abstract: Energy levels of a two-body atomic system in an external homogeneous magnetic 

field can be presented in terms of magnetic moments of their components, however, those 
CO , magnetic moments being related to bound particles differ from their free values. Study of 

bound g factors in simple atomic systems are now of interest because of a recent progress in 
experiments on medium Z ions and of a new generation of muonium experiments possible 
with upcoming intensive muon sources. We consider bound corrections to the g factors in 
^ , several atomic systems, experimental data for which are available in literature: hydrogen, 

■ helium-3 ion, muonium, hydrogen-like ions with spinless nuclei with medium Z. 
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, [Traduit par la redaction] 

Ph! 1. Introduction 

High-resolution spectroscopy of hydrogen and hydrogen-like atoms provides significant amount of 
(~| , data for precision QED tests and determination of basic fundamental physical constants such as the 

Rydberg constant Roc, the fine structure constant a, etc. However, not only the spectrum of a free 
hydrogen-like atom can be studied, but also its interaction with an external field and in particular with 
a homogeneous magnetic field. From experimental and practical point of view two kinds of two-body 
systems are of particular interest. Those are muonium and an hydrogen-like ion with a spinless nucleus 
at not too high Z. 

g factor of a bound lepton (electron or muon) can be presented in the form [1] 



g = 2-{l + a + b), (1) 

where 1 stands for the Dirac value, a is anomalous magnetic moment of a free lepton and an addi- 
tional correction b is due to the binding effects. The problems for calculation of b in muonium and an 
hydrogen-fike ion are similar but they have some different significant features. 
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2. Motivation: muonium 

Properties of the hydrogen atom are affected by effects of the proton structure which limit accuracy 
of any theoretical prediction. The muonium atom is, in contrast to hydrogen, a pure leptonic system 
suited for an accurate QED test. A value of its hyperfine structure has been measured at LAMPF with 
a great accuracy [2] 

i/hf,(exp) = 4 463 302.776(51) kHz . (2) 

The experimental uncertainty of Vhis is only 1.1 • 10 ~^ in fractional units. The calculation has not been 
able to reach this level yet. The theoretical problems are: 

• It is not possible to do any exact calculations. One have to expand either in a (i.e. to select 
diagrams with a different number of QED loops), or Za (i.e. over the binding effects) or m/M 
(i.e. over recoil effects). A possible contribution of presentely uncalculated terms of higher order 
has been estimated as 5 • 10~^ [3, 4]. 

• One has to take into account hadronic effects (see e.g. [3, 5]). 

• It appears that a largest problem of the theoretical calculation is an experimental uncertainty. 
The leading contribution into the ground state hyperfine splitting (so-called Fermi energy) 



Vp = -ir{Za) Z'-R^c f — 



(1 + flu) (3) 



contains the electron-to-muon mass ratio me/m^. The most accurate determination of the mass 
ratio has been performed at LAMPF [2] studying Breit-Rabi levels of the muonium ground state 
in presence of homogenous magnetic field. For any successful interpretation of the experiment, 
g factors of muon and electron in muonium are needed. 

Special features of the theory of g factors of electron and muon in muonium are: 

• It is necessary to find a theoretical expression for the bound g factors with an uncertainty below 
0.01 ppm. 

• All small parameters, a, Za and me/m^, are of about the same value (a = Za ^ 1/137 and 
me/ and hence QED, binding and recoil effects are equally important. 

The theory up to the third order in either of these three small parameters has been known up to now 
and our target is a complete evaluation of fourth-order corrections. 

3. Motivation: Hydrogen-like ions witli spinless nuclei 

Recent Mainz-GSI experiment [6, 7, 8] have provided precision measurement of 

• ion cyclotron frequency 

UJion = ^^TT^B (4) 

• and spin precession frequency 



???? NRC Canada 



Karshenboim and Ivanov 



3 



with uncertainty below 1 ppb. A comparison of these two frequencies leaves us with two values: of the 
g factor of a bound electron and of the electron-to-ion mass ratio. Accuracy of the latter is essentially 
the same as of determination of the electron-to-proton mass ratio in Ref . [9] and in fact the uncertainty 
for nie/Mion derived from Ref. |9| is bigger than both the experimental uncertainty for frequencies 
and the theoretical uncertainty for the bound g factor (see e.g. [10]). The comparison of experiment 
with theory provides the most accurate determination of the electron-to-proton mass ratio and a test 
of QED calculations. The latter is possible if the experiment on two different ions is performed and 
hence a ratio of their g factors is free of any uncertainty due to the proton-to-electron mass ratio. 
Published experimental data are related to the ion of carbon- 12 [6], recent efforts has been directed to 
the oxygen-16ion [8]. Precision measurements with other ions are also possible and expected. 

The theory must target a level of accuracy as good as few parts in 10^". The details of theory differ 
from muonium: the recoil effects are pretty small and less important than the QED contributions. The 
most important problem is to take into account the binding effects. The latter are less important for 
lighter H-like ions, namely for hydrogen- like ions of helium-4 and beryllium- 10. 



4. Muonium: calculations 



The "old" theoretical expressions for the ground state in muonium (up to the third order) are of the 
form (see e.g. [11]) 



5; (up to 3rd) = gf^ 
and 

5;(upto3rd) = 



1 



{Za) 



1 - 
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+ 
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(6) 



(7) 



Although Z = 1 for muonium we keep it to simplify classification of contributions. 
We present here a complete result for the fourth order corrections 



5; (4th) 



■,(0) 



{Zaf 
12 



iZaf{l + Z) 



- (0.289...) 



irif 
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(8) 



and 



<(4th) = g^^ 



= «(0) . } - 



a{Za){l + Z) 



TO,, 



+ 



a'^{Za) TOe 97 



6tt 



108 
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(9) 



Most of terms in Eqs. (8), (9) are of the kinematic origin and we have systematically checked them. All 
of those kinematic terms can be derived by several methods applied originally for derivation of the third 
order corrections which are all of the kinematic origin (cf. [11]). ^on-kinematic corrections include a 
high-order Breit contribution to g'^ [12] and a screening contribution to g'^j^. The latter was calculated 
using of the corrections to the Dirac wave function due to hyperfine interaction [13] and due to external 
magnetic field [14]. The result agrees with calculations for the screening tensor in Refs. [15, 16]. 

The results of a re-evaluation of LAMPF data [17, 2] taking into account Eqs. (8), (9) is shown in 
Table 1. 
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Table 1. Re-evaluation of LAMPF data. The fractional corrections Je,^ to the derived value of the ratio of the 
muon magnetic moment and proton magnetic moment are related to the fourth order contribution to g'e.M found in 
this work. 



Experiment 


LAMPF, 1982, [17] 


LAMPF, 1999, [2] 




3.183 290 0(13) 


3.183 289 03(38) 




3.183 3461(13) 


3.183 345 14(38) 


Se 


-0.4 • 10"** 


-0.3 ■ 10-" 


5, 


1.2 • 10"* 


1.2 • 10"^ 



5. H-like ions with spinless nuclei: calculations 

We present a theoretical expression for the binding correction to the electron g factor in the form 

b = 6kin + brel + &rec + ^rad + Kucl ■ (10) 

The kinematic term was discussed above in the case of muonium. The relativistic corrections are known 
exactly thanks to Breit [12]. The nuclear corrections were considered in [18, 1, 19]. For Z < 10 an 
analytic expression [1] 

^uci = ^ (^a)^ (miijv)' (11) 

has an adequate accuracy. 

The recoil effects are small enough. We have performed an eviduation based on the Grotch equation 
and confirm here the leading term 

1 771 

6. ec = -^(^ar^, (12) 

obtained recently in Ref. [20]. There were also numerical calculations by Shabaev and Yerokhin [22] 

and we apply their data to obtain our final results. 

The essential binding effects contribute into the radiative corrections 

&rad = byp + h\MK + &SE + ^21oop ■ (13) 

The vacuum polarization effects were studied in Ref. [18] numerically and in Ref. [14] analytically. 
The Wichman-Kroll contribution was investigated in Ref. [1]. The other light-by-light contributions 
are needed more study [21]. The two-loop contributions have not been yet calculated and only some 
conservative estimations are possible. The biggest contribution and the biggest uncertainty in &rad 
comes from the one-loop electron self energy. Fractional uncertainty of numerical data for 6se(^) 
encreases at lower Z and one can impove accuracy of low Z results after properly fitting data. We fit 
here numerical data from Ref. [18]. An essential improvement for Z = 2 and Z = 4 has been achieved. 
After our analysis has been completed a paper [23] with data more accurate than thoes in Ref. [18]. 
Results of our fitting are in fair agreement with new numerical calculations. 

The final results for ions of interest are collected in Table 2 below. Note, that the theoretical accu- 
racy for the hydrogen-Uke ions of helium-4 and berylhum-10 is higher than for carbon and oxygen and 
from theoretical point of view they are better suited for determination of the electron-to-proton mass 
ratio. 

A comparison of carbon-to-oxygen and oxygen-to-oxygen leads to 

5(^2c5+)/^(i6o7+) = 1.0004972734(9) , (14) 



???? NRC Canada 



Karshenboim and Ivanov 



5 



Table 2. The binding contributions to g factors of some hydrogen-like ions. 



Ion 



b 

[10-8] 




-70 948.84(3) 
-283 865.0(2) 
-638 857.4(5) 
1 136 142.5(8) 
1 136 141.9(8) 



1607+ 
1807+ 



and 



O7+)_g(i8O^+) = -1.2-10 



9 



(15) 



Equation (14) can be used to test QED calculations, while Eq. (15) is rather useful to check for experi- 
mental systematical errors. The experimental result for the ratio in Eq. (14) is 



where we take into account a possible correlation between systematic errors for the carbon experiment 
[6] and the oxygen measurement [8]. The experimental value is in a fair agreement with the theoretical 
prediction (14). 

6. Muon-to-electron mass ratio 

Several results of determination of the muon-to-electron mass ratio are collected in Table 3. The 
improvement of theory for the g factors of a bound electron and muon in the muonium atom is not that 
important for present level of accuracy. However, a number of promising projects for intensive muon 
sources for experiments in particle physics are now under consideration [36]. The accurate theory of 
the g factors in muonium will be useful for upcoming experiments. 

7. Proton-to-electron mass ratio 

Determination of the proton-to-electron mass ratio is summarized in Table 4 below. Presently the 
most accurate value of the electron-to-proton mass ratio comes from an experiment related to the 
g factor of a bound electron. Theory and experiment contribute into uncertainty at the same level. 
Experiments with lighter ions might reduce the theoretical uncertainty substantially. 

8. Tlie fine structure constant 

The mass ratios of electron-to-muon and electron-to-proton are considered above. Both are im- 
portant for precision determination of the fine structure constant a. A precision value of the latter is 
strongly needed for any accurate QED test. 

Comparison of theory and experiment for the muonium hyperfine structure interval allows to de- 
termine a if a value of the muon mass is known in a proper units, e.g in unit of the electron mass, or 
the muon magnetic moment in unit of the magnetic moment of electron or proton. 

A promising method to determine the fine structure constant has been developing by Chu and 
coworkers [45]. It is based on the photon recoil spectroscopy and offers a precision value of /i/Matom- 
The experiment in Ref. [45] deals with cesium atoms, but some other atoms are also used in several 



g(i2C''^+)/,g(if'0^+) = 1.000 497 2731(15), 



(16) 
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Table 3. Muon-to-electron mass ratio. The results marked with f differ from originally published. The result of 
Beltrami et al. is corrected accordingly to Ref. [4]. The result of Casperson et al. is presented accordingly to 
Ref. [24]. The results marked with f are found with taking into account fourth order corrections to the bound g 
factors of electron and muon. The CODATA result is based on some result derived from Is HFS Mu with some 
overoptimistic estimation of theoretical uncertainty. The result from Is HFS Mu quoted here is calculated with 
estimation of theoretical uncertainty as explained in Ref. [3] and with use of 1/a = 137.035 999 58(52) from 
anomalous magnetic moment of electron [25, 26]. 



Method 






Ref. 


Precession of ii^ in water 


206, 


.768 60(29) 


Crow ct al., 1972, [27] 


Precession of iJ.^ in Br2 


206, 


.768 35(11) 


Klempt et al., 1982, [24] 


(5-2) of/x- 


206, 


.771 4(21) 


Bailey et al., 1977, [28] 


(5-2) of/x+ 


206, 


.7668(20) 


Bailey et al., 1977, [29] 


(g-2) ofp+ 


206, 


.76870(27) 


Brown etal., 2001, [30,31] 


'i<k,/2 — 2p3/2 transitions 








in muonic ^''Mg and ^**Si 


206, 


.767 94(64)* 


Beltrami et al., 1986, [32] 


Breit-Rabi for Mu 


206, 


.768 18(54)* 


Casperson et al., 1977, [33] 


Breit-Rabi for Mu 


206, 


,768 22(8)+ 


Mariam et al., 1982, [17] 


Breit-Rabi for Mu 


206, 


,768 283(25)t 


Liu et al., 1999, [2] 


Is HFS Mu 


206, 


,768 283(11)* 


Liu et al., 1999, [2] 


Is — 2s in Mu 


206, 


.768 38(17) 


Meyer et al., 2000, [34] 



Adjustment 206.768 265 7(63) CODATA, 1998, [35] 



Table 4. Proton-to-electron mass ratio. The result for '"'Be"'" differs from the original value in Ref. [37] because 
we have applied a more accurate theory from Ref. [38] and a more accurate value of the beryllium-9 mass [39]. 
We present two theoretical evaluation of experiment [6]. The value marked with f has been obtained in this work. 



Method 




nip /nie. 


Ref. 


g factor in ^Be"*" 


1836. 


,152 92(29)* 


Wineland et al., 1983, [37] 


5factorini2c5+ 


1836. 


152 673 3(14) 


Haffner et al., 2000, [6] 








Beier et al., 2002, [40] 


g factor in ^'^C^+ 


1836. 


152 673 6(13)t 


Haffner et al., 2000, [6] 


Precession of e and p 


1836. 


152 667 0(39) 


Farnham et al., 1995, [9] 


Precession of e and 


1836. 


152 680(88) 


Gabrielse et al., 1990, [41] 


Gross structure in H and D 


1836. 


152 667(85) 


de Beauvoir et al., 2000, [42], 








Huber et al., 1998, [43] 


Antiprotonic He 


1836. 


152 58(24) 


Hori etal., 2001, [44] 


Adjustment 


1836.152 667 5(39) 


CODATA, 1998, [35] 
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experiments by other groups, order to obtain precise value for a, one can combine h/M^^om with a 
precision value of the Rydberg constant 

- T ~r ' ^^^^ 

the cesium-to-proton mass [46] and the electron-to-proton mass ratio. 

The presently most accurate value of the fine structure constants comes from study of anomalous 
magnetic moment of electron [25]. It used to be quoted as a QED value. However, the major uncertainty 
of this value of a comes from not perfect understanding of a motion of a light particle (i.e. electron) in 
a Penning trap [26]. It is important that the value of the electron-to-proton mass ratio from the bound 
g factor is free of problem of an electron in Penning trap in contrast to the former best value of the ratio 
from University of Washington [26]. That makes a value of a from Chu experiment really independent 
from those from g — 2oi electron. The recent preliminary result by Chu has a relative uncertainty as 
low as 7 • 10~^ and some progress is still possible [45]. 
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